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I. INTRODUCTION 

Investigation of the critical behavior of mag- 
netic systems is usually made harder, compared to 
the translationally-invariant case, when one considers 
quenched disorder. This is because one gets inherent, 
randomness-induced, fluctuations in experimentally- or 
numerically measurable quantities, whose effects are of- 
ten difficult to separate from those of purely thermal ori- 
gin, or connected with finite-size scaling and crossover 
phenomena. Therefore, redoubled interest is attracted 
when exact results (either rigorously proved or conjec- 
tured) are put forward in the context of disordered spin 
systems. Suitable tests can then be devised in which, 
by taking advantage of the proposed exact relationships, 
one attempts to provide a clearer physical picture of the 
problem at hand. 

Here we study two-dimensional ± J Ising spin glasses, 
i.e. Ising spins interacting via nearest-neighbor bonds 
of the same strength and random sign. While no spin- 
glass ordering arises at T ^ for equal concentrations 
of ferro- (p) and antiferromagnetic (1 — p) bonds, in the 
asymmetric case p ^ 1/2 one can have long-range order 
for suitably low concentrations of frustrated plaquettes. 
Accordingly, for p not far from unity a critical line on the 
T — p plane separates paramagnetic and ferromagnetic 
phases. 

A number of exact results have been derived along a 
special line in the T—p plane, the Nishimori linei^; also, 
the exact location of a multicritical point, the Nishimori 
point, believed to be at the intersection of the critical 
boundary with the Nishimori line, has been predicted^. 

In this paper, we investigate the multicritical behavior 
at the Nishimori point. We use numerical transfer-matrix 



methods to calculate probability distributions and asso- 
ciated moments of spin-spin correlation functions, and 
try to ascertain whether their properties may reflect con- 
straints imposed by conformal invariance requirements. 

In section |n] we recall some exact statements and nu- 
merical results, pertaining to the Nishimori line and the 
location of the Nishimori point. In section IIIII we de- 
scribe the numerical techniques to be used, and provide 
a test of their correctness by performing calculations on 
special unfrustrated systems. In section ITvl we display 
and analyse results from numerical calculations of corre- 
lation functions on strips, concentrating on the shape and 
properties of the corresponding distributions, especially 
as regards: (i) correlation-function equalities stemming 
from the gauge symmetry obeyed by Ising spin glasses, 
and (ii) conformal-invariance requirements. Section IVl is 
dedicated to an approximate scaling scheme which pre- 
serves some of the essential symmetries obeyed along the 
Nishimori line. Finally, in section lVll concluding remarks 
are made. 



II. ±J ISING SPIN GLASSES 

We consider Ising spin-1/2 variables on sites of a 
square lattice, with couplings between nearest-neighbors 
i and j drawn from the binary distribution: 

P(Jy) = P 5(J ij - J ) + (1 -p) S{Ja + Jo) • (1) 

The Nishimori line (NL) is defined by the following rela- 
tionship between temperature T and positive-bond con- 
centration p: 

e -2j /T _ y (Nishimori line, p > -) . (2) 
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Along this line, the configurationally-averaged internal 
energy is analytic and can be calculated in closed form 1 . 
Also, a number of exact relationships holds between mo- 
ments of correlation function distributions, as well as be- 
tween order parameters and their derivatives. Of particu- 
lar interest here will be the following property, which has 
been shown to hold on the NL, for correlation functions 
between Ising spins Oi 
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where angled brackets denote the usual thermal average, 
square brackets stand for configurational averages over 
disorder, and I = 0, 1, 2, ... . 

It is widely believed that, within the ordered phase, the 
NL separates a high-temperature region, dominated es- 
sentially by pure-system behavior, from one where zero- 
T effects play the leading role. A doubly-unstable (mul- 
ticritical) fixed point is then expected to exist on the 
ferro-paramagnetic phase boundary. It was argued^ that 
the location of this, also known as Nishimori point (NP), 
should coincide with the intersection between the criti- 
cal boundary and the NL. We shall take this view in the 
following. 

Recently it was predicted^! that, on a square lattice, 
the NP should belong to a subspace of the T — p plane 
which is invariant under certain duality transformations. 
For ±J Ising systems, the invariant subspace is given 
byM: 

p log(l+e- 2Jo/T ) + (l-?j) log(l+e 2J °/ T ) = - log 2 . (4) 

The NP is thus predicted to be at the intersection of 
Eqs. and JU, namely p = 0.889972 T/J = 
0.956729 • • • . This agrees well (though, in some cases, 
it is slightly outside estimated error bars) with avail- 
able numerical results, from which p is given respec- 
tively as: 0.886(3) (series)^, 0.894(2) (zero-T calcu- 
lations, assuming the phase boundary to fall verti- 
cally from the NP)2>; ~ 0.885 (exact combinatorial 
work)itt; 0.8905(5) (transfer-matrix scaling of correla- 
tion lengths)ii; 0.8906(2) (transfer-matrix scaling of 
domain-wall energies)^; 0.8907(2) (mapping into a net- 
work model for disordered non-interacting fermions, via 
transfer-matrix} 13 ; 0.8894(9) (Monte-Carlo analysis of 
non-equilibrium relaxation)^. 

Though the pure-system critical point also belongs to 
the subspace given by Eq. @J, it is known£ that inter- 
preting that equation as the exact form of the critical 
boundary (at least above the NP) brings problems of its 
own. Indeed, at p = 1, the reduced slope of the phase 
diagram is predicted^ to be exactly: 
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While numerical transfer-matrix calculations are in good 
agreement with this, giving respectively 3.25(ll)ii and 



3.23(3)i£, Eq. ® yields 2 + y/2 = 3.4142 • • ■ , 6.4% in 
excess of Eq. Such discrepancy is to be compared 
to the disagreement between the above prediction for the 
NP and the various central estimates quoted there, which 
never exceeds 0.5% on either side. 

Thus, although the conjectured location of the NP (to 
be referred to as CNP) given by Eqs. |J2Jl and Q may turn 
out not to be exact, it is certainly a good approximation, 
and will be taken as the starting point in what follows. 

On the other hand, conformal invariance properties^ 
are known to hold at the critical point of pure two- 
dimensional magnets, and evidence has been provided 
that they are present in random systems (at least, unfrus- 
trated ones) as well, namely: random-bond IsingiiA&i^, 
random-bond g-state Potts models 2 ^, and random trans- 
verse Ising chains at T = (equivalent to the two- 
dimensional McCoy- Wu model)2i. As regards the two- 
dimensional Ising spin glass, early attempts to ap- 
ply conformal invariance ideas at the NPii were un- 
dertaken in the context of testing whether the ferro- 
paramagnetic transition was in the same universality 
class as random percolation, as suggested by series work 8 . 
While extrapolation of the uniform zero-field suscepti- 
bility gave the exponent ratio 7/1/ = 1.80(2), consis- 
tent with (j/v) P = 43/24 = 1.7917- •• of percolation, 
the exponent-amplitude relation 2 ^ L/tt^l = V (£l is the 
correlation length on a strip of width L, at criticality) 
gave rj = 0.182(5), thus excluding the percolation value 
Vp = 5/24 = 0.208333 ••• . However, it must be noted 
that the calculated values of 7/1/ and 77 are not mutually 
excludent, via the scaling relation 7/1/ = 2 — 77. This 
indicates that a reanalysis of the validity of conformal 
invariance at the NP is in order (though the connection 
to percolation can probably be ruled out, as shown by 
later numerical dat a 12 ' 13 ). Also, one can devise more 
stringent tests of conformal invariance requirements than 
that given by the exponent-amplitude relation. Indeed, 
the correlation length entering that relation is the single 
parameter governing the asymptotic (exponential) decay 
of correlations, whereas for short distances conformal in- 
variance implies specific functional relationships^. This 
fact has been exploited in Ref. 0, where assorted mo- 
ments of the correlation function distributions were fitted 
to a form suggested by conformal invariance arguments, 
thus enabling the extraction of the corresponding expo- 
nents from relatively short-range correlation data. 

Our goals here are: (1) to verify the extent to which 
the consequences of Eq. <EJ| are directly reflected in the 
probability distributions of the Cy; and (2) to probe the 
angular dependence of correlation functions, in order to 
test whether they obey the predictions of conformal in- 
variance. 



III. CORRELATION FUNCTIONS ON STRIPS 

We apply numerical transfer-matrix (TM) methods to 
the spin-1/2 Ising spin glass, on strips of a square lattice 
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of width 4 < L < 12 sites. We have calculated correlation 
functions C xy between spins separated by x lattice spac- 
ings along the strip, and y in the transverse direction 1 ^ 
such that R = (x 2 + y 2 ) 1 ^ 2 is typically of order L. On 
account of periodic boundary conditions across the strip, 
we need only consider < y < L/2. By iterating the TM 
on sufficiently long strips, we have accumulated enough 
non-overlapping samples of correlation functions (usually 
TV = 10 5 — 10 7 ), in order to produce clean histograms, 
P(C xy ), of occurrence of C xy . We have used a linear bin- 
ning, dividing the [—1,1] interval of variation of C xy into 
N hin = 10 3 equal bins of width S = 2 x 1CT 3 . Further- 
more, we have used a canonical ensemble (instead of the 
usual grand-canonical scheme) for the bond distribution, 
so that overall concentration fluctuations are kept to a 
minimum. 

Correlation functions are calculated along the lines of 
Section 1.4 of Ref . E^L with standard adaptations for an 
inhomogeneous system^. For two spins in, say, row 1, 
separated by a distance R along the strip, and for a given 
configuration C of bonds, one has: 

{VqVr/C - Z 7 R1 x i 

Eaoafl^K) (IIi=0 T i) ^r) 

(6) 

where <jq = {<Tq . . . <Jq } an d correspondingly for ctr; the 
bonds that make the transfer matrices % belong to C. For 
pure systems the 2 L — component vectors ip, ip are deter- 
mined by the boundary conditions along the strip; for 
example, the choice of dominant left and right eigenvec- 
tors gives the correlation function in an infinite system 2 ^. 
Here, one need only be concerned with avoiding start-up 
effects, since there is no convergence of iterated vectors. 
This is done by discarding the first few hundred iterates 
of the initial vectors, Vq , vo. From then on, one can 
shift the dummy origin of Eq. U3 along the strip, taking 
-0 (ip) to be the left- (right-) iterate of (v ) up to 
the shifted origin, and accumulating the corresponding 
results. In order to avoid spurious correlations between 
the dynamical variables, the respective iterations of ip 
and ip must use distinct realizations of the bond distri- 
bution. 

The correctness of the procedures just described can 
be checked by testing them on the ferromagnetic (FM) 
random-bond case with a distribution of couplings given 

by 

P(Jij) = \Wij-Jo)+${Jii-rJ )), < r < 1 (FM). 

(7) 

In this case, the exact critical temperature /3 C — l/fcsT c 
is know n 24 i 25 to be 

sinh(2/3 c Jo) sinh(2/3 c r Jo) = 1 (FM). (8) 

It is known from Monte-Carlo work^Sl on L x L lattices, 
L < 1024, that the critical correlation functions of the 
FM random-bond model are numerically very close to 
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Figure 1: Triangles: averaged correlation functions for FM 
model at criticality (see Eqs. iQand (01 ). Error bars are twice 
the standard deviation among three independent runs. Con- 
tinuous line connects pure-system results for infinite strips, 
i.e. calculated with ip, ij) of Eq. as the dominant left and 
right eigenvectors of the TM. 



those of a pure Ising system at its own criticality; for 
R/L < 1 the discrepancy is always under 4%, and ap- 
proaches zero as R/L — > 0. 

We considered a strongly disordered system with r — 
1/4, and calculated the averaged correlation functions 
[(co0\r)] (first moment of the distribution) on strips of 
width 4 < L < 16. For each strip width we made three 
independent runs, each long enough to collect N non- 
overlapping samples of correlation functions for a; = L/2, 
y = (N = 10 s for L < 12, 5 x 10 4 and 4 x 10 4 respec- 
tively for L = 14 and 16). The results are shown in 
Fig. where error bars are twice the standard devia- 
tion among runs^. One can see that the near identity 
with pure-system critical correlations is recovered (for all 
points in the Figure, central estimates differ from their 
pure-system counterparts by less than 1%) . We conclude 
that our calculational methods, in particular the choice 
of vectors ip and ip, are suitable for the description of 
correlation functions in disordered systems on strips. 



IV. NUMERICAL RESULTS FOR ±J ISING 
SPIN GLASS ON STRIPS 

We now return to spin glasses. We first note that the 
pairing of successive odd and even moments predicted in 
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Eq. implies that 
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(l-Cy)P(Cy) dG« = .(9) 



Since this holds for any £, P'(C xy ) = (1 - C xy ) P(C X?/ ) 
must be an even function of C xy , everywhere on the 
NL. For the time being, we shall only consider y — 
0. Fig. El where results at the CNP for L = 4, 8, 
and 12, (x,y) = (L/2,0) are displayed, illustrates that 
such parity property is present, apart from small fluc- 
tuations caused by finite-sample and binning effects. 
The quantity A = A(C) = P'{C xy ) - P'{-C xy ) pro- 
vides a quantitative measure of these deviations. With 
A^ = [(N hin /2)- 1 Z 1 c=a [A(C)] k } 1 /\ we find that: (i) 
for N = 2 x 10 6 , IA 1 ! < 1CT 6 and A 2 ~2x 1CT 5 , con- 
stant against L; and (ii) for fixed L = 4 (where lattice 
structure has more pronounced effects, see Fig. thus 
one would expect the most unfavorable environment for 
measuring such small fluctuations) and 10 4 < N < 10 7 , 
A 2 varies roughly as N~ x / 2 . The latter is to be expected 
if deviations from parity arise from (randomly) incom- 
plete sampling. Therefore, we can be confident that the 
property given in Eq. Q is satisfied by our data, except 
for small random deviations whose origin is well under- 
stood. For N = 2 x 10 6 samples and 4 < L < 12 the 
quantitative effect of such deviations on the moments of 
the distribution is that, for < I < 3, Eq. is satisfied 
to within one part in 10 4 . This holds not only for the 
y = data shown above, but for all relative positions 
(x,y) investigated here (to be discussed below). 

As regards conformal invariance, we first recall that, 
for pure Ising systems on a strip of width L with periodic 
boundary conditions across, the following result holds at 
criticalityifi: 



xy 



n/L 



(sinh 2 (7ra/L) + sin 2 (vry/L)) 



1/2 



, ?7=l/4, 
(10) 

where the proportionality factor can be obtained from 
the exact square-lattice (L, R — » oo, R -C L) resultSl, 
C R = 0.703 38/i? 1 / 4 . Though strictly speaking Eq. QDJ 
is an asymptotic form, we have checked that discrepancies 
are already very small at short distances: for an L = 
20 strip, the largest deviation to be found between that 
and numerically-calculated values is 1.3% for (x, y) — 
(1,1), while it remains below 0.5% for all other cases, 
approaching zero very fast: for (x,y) = (3,3) it is 0.2%. 
The picture is the same already for smaller strip widths, 
e.g. L — 10. 

Turning to disordered systems, a version of Eq. 
was used at the NP in the strip calculations of Ref. | 
Moments of order i = 1, 2, • • • 8 of the correlation function 
distribution, for y — and varying x < L, were numeri- 
cally calculated and fitted to the form Eq. (1101 . with the 
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Figure 2: Data taken at the CNP. P'(C) = (1 - C)P(C), 
where P(C) is the normalized histogram of occurrence of cor- 
relation function C. Strip widths L = 4, 8 and 12 (the latter 
two vertically shifted respectively by 0.002 and 0.004, to avoid 
superposition); spin-spin relative positions (x = L/2,y = 0). 
N = 2 x 10 6 non-overlapping samples, for all cases. 
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Figure 3: Differences 2?, between symmetrized distribu- 
tion functions P'(C xy ) - P'{C x , y ,) ( P'(C) = (1 - C)P(C), 
where P(C) is the normalized histogram of occurrence of 
correlation function C). V x = P'(C 15 ) - P'(C 23 ), V 2 = 
P (C15) — P (C&i); insert(same scale as main plot): T>3 = 
P (C23) - P'(C 3 2). Data taken at the CNP for strip width 
L = 10, N = 2 x 10 6 non-overlapping samples in all cases. 
Dashed horizontal lines at ±A 2 (rms deviation of P'(C) from 
parity for N as given above, see text). 
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corresponding iji to be extracted from the fitting proce- 
dure. For i — 1 their result is T]i — 0.1854 (error bar 
estimated as ~ 1% 12 >), which agrees with, and is more 
accurate than, that obtained from the asymptotic decay 
of correlations, namely 0. 182(5)^. 

Here we shall try to capture variations in the correla- 
tion function distribution and its moments, against vari- 
ations in the argument of Eq. ijlOj) . z = (sinh 2 (7nr/L) + 
sin 2 (7rj//L)) 1 / 2 . As (i) the coordinates (x,y) only take 
discrete values, and (ii) we have no guide as to the specific 
functional dependence of P(C xy ) on z, we shall at first 
compare distributions corresponding to points for which 
the argument is nearly equal, in contrast to those corre- 
sponding e.g. to the same distance, but with significantly 
different arguments. We take L = 10 and the points 
(x, y) = (1,5), (2,3), (3,2) and (5,1), for which one has 
W*16 = 1-001; z 51 /z 15 = 2.1922; z 32 /z 23 = 1.1773. 

In Fig. we show the differences between pairs of sym- 
metrized distributions P'(C xy ): V x = P' (C 15 ) - P' (C 23 ) , 
V 2 = P'(C 15 ) - P'(C 51 ), V 3 = P'(C 23 ) - P'(C 32 ). In- 
deed, the differences between distributions behave quali- 
tatively as one would expect, should the dependence on 
x and y be only through z. The horizontal lines at ±A 2 
(rms deviation of P'(C) from parity for N = 2 x 10 6 as 
is the case here) show that, for T>\ the 0.1% difference 
in the respective arguments of z gives rises to effects of 
the same order of magnitude as those arising from each 
individual distribution's deviations from parity. 

A more quantitative perspective can be obtained by 
analysing the dependence of assorted moments rrii of the 
distribution against z. In Fig. we show odd moments 
of the correlation function distribution at the CNP, for 
L = 10 and 1 < £ < 3, < j/ < 5. These correspond to a 
region where the argument z is strongly influenced by y, 
thus it is especially convenient for the discussion of the 
angular dependence implied in Eq. lll(H . 

From the data displayed above, the overall conclu- 
sion is that conformal invariance holds at the NP, in the 
strict angle-dependent form corresponding to Eq. ifTfil) . 
Note that the point (x,y) = (1,1), z = 1.0943, for 
which the discrepancy between numerically-calculated 
value and asymptotic expression of correlation function 
is the largest for pure systems, also corresponds to the 
largest deviation from least-squares fitting lines in the 
disordered case. 

We now turn to the question of whether the location 
of the NP, as conjectured in Refs.00, may be exact. We 
investigate the numerical values of the exponents Tj2j+i, 
as obtained from the moments of the correlation-function 
distribution. We have performed least-squares fits of data 
for fixed L = 10 with the grid of z— values shown in 
Fig. 0] (henceforth called z— fits), to the form suggested 
by Eq. lUf))) for this case, namely rri2j+i ~ z~ V2i+1 . We 
have scanned the NL in the immediate neighborhood of 
the CNP, 0.8886 < p < 0.8914, so as to include the un- 
certainty ranges of some recent estimates for the location 
of the N P 12 i 13 i 14 . We found that the quality of fits, as 
given by the corresponding chi-square, does not change 
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Figure 4: Double-logarithmic plot of odd moments of 
the correlation-function distribution P(C xy ) against z = 
(sinh 2 (7ra;/L) + sin 2 (ny/L)) 1/2 ; 1 < x < 3, < y < 5. 
Straight lines are unweighted least-squares fits to data. Data 
taken at the CNP for strip width L = 10, N = 3 x 10 6 non- 
overlapping samples in all cases. 



Table I: Estimates of exponents r/2j+i, from least-squares 
fits of averaged odd moments of correlation-function distribu- 
tions; z— fit: data for L — 10 and 1<£<3, 0<y<5, 
assuming m2j+i ~ z~ V2 ^ +1 - L = 4 — 12: data for x = L/2, 
y = 0, assuming rri2j+i ~ L~ V2j+1 . First and second 
columns: data taken at the CNP. Third and fourth columns: 
pure system at criticality. Uncertainties in last quoted digits 
are shown in parentheses. 
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noticeably along that interval. Our estimates are dis- 
played in Fig. where the pairs of dashed lines in the 
central region of the Fig ure show ranges of exponent es- 
timates given in Ref . ll2L namely T)2j+i = 0.1854, 0.2561, 
0.3015, 0.3354 for j = 0, 1, 2, 3 (error bars assumed to 
be l%i£). One sees that the results of z— fits are closest 
to those of Ref. Il2l at the CNP, not at the location of the 
NP predicted in that reference. 

In order to infer how much systematic error is implicit 
in the z— fitting procedure, we (i) performed least-squares 
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Figure 5: Results of least-squares fits of odd moments of the 
correlation-function distribution P{C xy ) to the form Eq. 1KH . 
with the r\i as parameters; L = 10, 1 < x < 3, < y < 5, 
N = 3 x 10 6 non-overlapping samples in all cases. Data taken 
along the NL, parametrized by ferromagnetic bond concen- 
tration p. Key to symbols is the same as in Fig. 0] (bottom 
to top: j = 0, 1, 2, 3). Error bars are of same order as sym- 
bol sizes. Vertical arrow: location of CNP. Heavy horizontal 
lines close to bottom show ranges of some recent estimates 
for location of NP. Long dashes: Ref. IT3 : full line: Ref. ITsI : 
short dashes: Ref. Q. Pairs of dashed lines show ranges of 
exponent estimates given in Ref. Il2l 



fits of data taken at the CNP for 4 < L < 12, x = L/2, 
y = against the form suggested by Eq. (II fl for that 
case, i.e. m2j+i ~ L~ r ' 2j+1 ; and (ii) took pure-system 
data at criticality, from the same grid of z— values used 
at the NP, and extracted the corresponding estimates of 
V2j+i from z— fits. The results of (i) and (ii) are shown 
respectively in the second and third columns of Table [I] 

From the analysis of pure-system data in Table Q] we 
conclude that, at least for j = and 1, the intrinsic 
inaccuracy of z— fits does not significantly compromise 
our final results. For these same values of j, we also find 
greater internal consistency between z— fit and L = 4— 12 
estimates at the CNP. This signals that finite-width ef- 
fects are essentially subsumed in the explicit L— depen- 
dence given in Eq (1101 . i.e. higher-order finite-size cor- 
rections presumably do not play a significant role. Of 
course, a decisive test of the latter statement would in- 
volve e.g. doing z— fits on significantly wider strips. 

For the moment, the data at hand allow us only to 
conclude that the conjectured exact location of the NP, as 
given in Refs.00, cannot be definitely ruled out. On the 
contrary, if we assume that the numerical values of the 
r\i given in Ref. XUh are the most reliable ones at present, 



then our results are in fact consistent with the conjecture. 



V. APPROXIMATE SCALING 

In order to further understand the properties (so far, 
numerically found) of correlation function distributions, 
in this section we introduce an approximate scaling trans- 
formation which incorporates the gauge symmetries un- 
derlying the NL . 

We recall, from real-space renormalization ideas, the 
concept of decimatior^^^: considering e.g. a square 
lattice, summing upon the degrees of freedom of spins 
on a given sublattice will produce a new square lattice 
with half as many sites as the original one, and a lat- 
tice parameter enlarged by a factor of y/2. In this con- 
text, the appropriate variable is the bond transmissivity 
t = tanh( J/T), which represents the spin-spin correla- 
tion function in a restricted (single-bond) ensemble. For 
a point on the NL, the constraint given in Eq. means 
that, with to = tanh(Jo/T), the probability distribution 
of t is 



1 



1 



V(t) = -(1 + 1 ) S(t - t Q ) + -(1 - t ) S(t + 1 ) . (11) 

Denoting by t i: i = 1,2,3,4 the transmissivities of 
bonds around a plaquette, the scaled transmissivity of 
the diagonal (effective) bond after decimation i a 2 ^^ 



t' 



1 + txtit^ti 



(12) 



We investigate how the distribution function of scaled 
transmissivities, V(t'), evolves under successive rescal- 
ings via Eq. l(T2ll . given that the original distribution 
obeys Eq. ifTTll . 

We first prove that the distributions generated by iter- 
ation of Eq. (11211 obey the NL symmetry given in Eq. 10, 
i.e. the NL is an invariant subspace of the transformation 
Eq. l(T2|l . The proof is by induction, as follows. 

(i) Direct examination of Eq. (1 1 XII shows that the mo- 



obey 



ments mT' of the original, i.e. zero-th order distribution, 



m 



(0) 
2^+1 



2£+l] (°) 



(0) 
l 2l+2 



2£+21 (°) 



(13) 



where the upper index (0) attached to square brackets 
denotes average over the distribution Eq. (1 1 111 . 



(ii) Assuming that 



rr> 



(n) 
21+1 



m 2£+2 ' whence 



[t 2e+1 (1 - t)] {n) = 0, where the upper index (n) denotes 
averages over the n-th order iterated distribution, one 
can show that 



£1*2 + ht4 



1 + ^1*2^3^4 



2e+l 



1 - 



£1*2 + £3*4 
1 + t\t2ty,t^ 



1 («) 



(14) 
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by expanding the RHS and invoking the statistical inde- 
pendence of the variables tj. 

In fact, the invariance property just proved is guaran- 
teed by Eq. 11211 . because the decimated spins are con- 
nected to an even number of bonds. Indeed, the gauge 
transformation originally introduced for the ±1/2 Ising 
spins, in the derivation of the NLi*fi>, could equally well 
have been applied to the transmissivities of the current 
approach, with the same final result. 

Analytical examination of the evolution of V(t) under 
iteration of Eq. l(T2|l shows the following relevant features: 

(1) Assuming that the distribution has a non-negligible 
weight near t — and making a self-consistent ansatz, in 
which V {n \t) ~ t a " , one obtains the fixed-point value 
a* = -1, i.e. V{t) ~ It]- 1 , t-> 0. 

(2) Assuming now that the distribution has a non- 
negligible weight near t = 1, and making an ansatz in 
which, with 5 = 1-t, Q (n) (*) = T {n \l - t) ~ 

one obtains the fixed-point value (3* = —1, i.e. V(t) ~ 

The analytic predictions for the exponent values a*, 
(3* were determined by the limiting forms of Eq. 1121) for 
small ti or small 1 — respectively. Of these, the first is 
the less reliable, since the simple transformation Eq. l(T2|l 
captures less well the key subtle effects needed there, of 
removal of order by frustration, than the simpler charac- 
teristics of the ferromagnetic state involved at small 1— ti. 
So we expect (and find, see below) that the predicted a* 
is much less reliable than the predicted (3*. 

In order to test the above predictions, one might work 
out the iterated distributions directly from Eq. ltl2ll . 
This, however, is feasible only for the first four iterations, 
since the number of distinct outcomes grows very rapidly 
from 2 of Eq. l(TT|l at order zero, respectively to 3, 5, 25, 
702 at subsequent orders (the latter figure, for the fourth- 
order distribution, reflects the binning of 25 4 = 390, 625 
results onto 10 3 bins of width 2 x 10~ 3 , while no bin- 
ning was used for the preceding iterations). The next 
step would involve considering of order 2 x 10 11 opera- 
tions, which might be reduced by perhaps one order of 
magnitude by considering symmetries. The alternative of 
considering a coarser binning structure (to be kept fixed 
for iterations of order > 4) was seen to introduce gross 
distortions of the symmetry reflected in Eq. (PJl, thus the 
presumptive advantages of exact recursion relations over 
the purely numerical techniques of section |^ would be 
lost. 

Before setting out to compare predictions (1) and (2) to 
numerically-collected data on strips, we mention that the 
soundness of the rescaling scheme can be roughly tested 
as follows. Following the usual renormalization-group 
ideas, one would expect the NP to be characterized by 
a scaling-invariant V*(t)^. In practice, one searches for 
equality between the j-th moments of the m-th and n-th 
order iterated distributions, expecting that the results (i) 
do not depend strongly on the choice of j, and (ii) should 
improve as m and/or n grows. Considering m = 0, one 
can parametrize the search for a fixed point distribution 



i — i — i — i i i i | 




Figure 6: Double-logarithmic plot of P(C) close to C = 0. 
Full symbols: C > 0; empty symbols: C < (divided by 1.10, 
to reduce superposition). Strip widths L = 8 (triangles)), 10 
(squares), and 12 (hexagons). For all cases, x = L/2, y = 0. 
Straight line corresponds to P(C) ~ | C | 1 / 3 . Data taken at 
the CNP. 




0.080.1 0.2 
1-C 



Figure 7: Double-logarithmic plot of P(C) close to C = 
1. Strip widths L = 8 (triangles)), 10 (squares), and 12 
(hexagons). For all cases, x = L/2, y = 0. Straight line 
corresponds to P(C) ~ (1 - C)~\ Data taken at the CNP. 
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through the value to in Eq. lllljl : for to = to, m j — m j - 
In this context, recalling the invariance of the NL under 
Eq. O and taking j = 1, we have p* = p(t$) = 0.9395 
and 0.9417, respectively for n — 3 and 4. An arf /ioc 
two-point extrapolation against 1/n gives a non-trivial 
limiting value, lim n ^ooP* = 0.959. This is to be com- 
pared with the presumed exact p — 0.889972 • • We 
conclude that the rescaling scheme is qualitatively cor- 
rect, though its quantitative predictions must be care- 
fully scrutinized. 

We have tested prediction (1) above against numerical 
data from strips. The results are displayed in Fig. El 
One can see that, although there are clear signs of a 
power-law divergence both below and above C = 0, the 
corresponding exponent is certainly not above 1/3. The 
slight asymmetry between data for C > and C < 
is entirely consistent with the general property that the 
symmetrized distribution P'(C) is an even function. In- 
deed, by plotting P'(C) we get two essentially identical 
branches. Turning now to prediction (2), the correspond- 
ing results are displayed in Fig. Now one can see that 
the form P(C) ~ (1 — C) _1 indeed captures the essential 
features of behavior. 



VI. CONCLUSIONS 

Our investigation has probed behavior at the Nishimori 
point (NP), where the infinite two-dimensional system is 
both critical and subject to the Nishimori gauge symme- 
try. Mounting, but incomplete, evidence has suggested 
that random critical two-dimensional systems share, in 
a statistical sense, the conformal invariance of their pure 
counterparts. So it was expected that at the NP, the cor- 
relation function distributions P(C) on strips would be 



the same at all (lattice) points (x, y) on any curve of con- 
stant | sinh(7r(a; + iy)/L) | (see Eq. itTTHl l. Evidence has 
been presented in section |^ that this strict form of sta- 
tistical conformal invariance indeed applies. This is espe- 
cially important here, since all numerical evidence avail- 
able so far, regarding conformal-invariance properties at 
the NP, points to values of e.g. the exponents rji and the 
central charge which are not obviously associated to any 
previously known universality classii*iSiH . Our data do 
not rule out the possibility that the recently-conjectured 
location of the NP, at the intersection of Eqs. @ and Q), 
is indeed exact. 

Among other results are the power-law divergencies of 
distributions P(C) near C = 1 and C = 0, which were 
first identified in the invariant distributions of the simple 
scaling theory (section^ and then confirmed by the strip 
scaling analysis (section^. Of the respective power laws 
(— 1, —1/3) only the first was correctly predicted because 
of the crudeness of the scaling transformation, Eq. Q2j| 
(which nevertheless maintains the gauge symmetry). 
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